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1. Introduction 

The needs of nuclear material engineering for intensive 

sources of nuclear radiation are growing all the time. Reactors 

are being constructed in which the materials and assemblies af 

various devices can withstand ever increasing fluxes af thermal 

and fast neutrons, which makes it passible ta study destruction 

of the fuel assemblies (Ref.l), obtain trans-uranium elemeuts 

(Ref, 2), investigate the properties of short-lived nuclei (Ref, 3). 

14 -2 -1 

Ten years ago a thermal neutron flux af about 10 cm sec was 

considered ta be very high whereas at present there are reactors 

16 —2 —1 

which produce short-time fluxes as high as ^ 10 cm sec • It 
may be supposed that the fluxes will grow t« still higher values. 

However, since each specific method of obtaining intensive 
nuclear radiations has its limitations these methods will change In 
quality. One af the new methods is the use of underground nuclear 
explosions for obtaining powerful neutron impulses (Ref .4). 

Studying the maximum potentialities of the used powerful 
neutron sources is a main theoretical problem. It will be shown 
below that for a given amount of energy which can be taken by one 
volume unit of the coolant there is a constant limit to which the 
neutron density in the core of the reactor can be increased, re- 
gardless of the known reactor type. 

The paper decels with pulsed reactors, reactors with fuel 
motion relative to the core, xoactors with transfer of heat from 
fuel to coolant and reactors with opposing coolant streams. 
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2, l\ilaocl roacboi'Q 

The pulsed roactox' koopa operating unitl it heats up to the li- 
mit imposed by its Q«iiS't£uotl*n and the materials . Then the reactor 
has to cool dovm before it again becomes operative (the heat loss 
during the operation is usually negligibly small). 

Apart from cooling the operation cycle of a pulsed reactor can 
be divided into three stages: build-up of power to the maximum level 
operation at this level, collapse of the power. As the reactor heats 
up during all three stages, the shorter the duration of each stage, 
the higher the maximum power. 

The rate of neutron density n growth (and, consequently, the 
rate of power growth) is limited by an exponential curve with a pe- 
riod. 

•T . _ ^ 

where (K«, -1) is the maximum possible supercriticality of the 
reactor (of infinite dimensions), V - velocity of fission-producing 
neutrons, £" -abaorbtion cross-section of these neutrons. With 
such neutron density growth when it reaches maximum the integral 
(in time) neutron density is equal to 

Jndt = (2) 

Accordingly, the heat accumulated per volume unit of the reactor 
will be equal to Q = ^f " 

fission energy, - neutron yield per fission event, - 

liberation of energy corresponding to the absorbtion of one neutron. 

The value of Q cannot exceed the permissible heat liberatioii 
per volume unit of the reactor The equality Q and is 

possible only with the zero duration of the operatlag duty and in- 
stanteneous collapse of power. Therefore the neutron density cannot 
exceed the value given T^y 

z ( 3 ) 

which means that in accordance with (1) 


.. ' (yf - 


>CU 
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(v) 

Let us represent as the number #f attme per volume unit of 

the reactor (N) multiplied by the nioximum energy which can be 
on the average imparted to each atom ( S ): 



If the construction of the reactor does not provide for the tran- 
sition of the core to the plasma state, the value of 6 must 
not exceed several tenths of an elecrron-volt • Assuming that 

we obtain 

14 

%ax 

The neutron flux depends on the energy of the neutrons which 
cause the bulk of fissions. In the thermal energy region the cross- 
section of many important reactions with neutrons follows appro- 
ximately the 1/V law, therefore "n” is in many cases a sufficient 
characteristic for thermal-neutron reactors. It is possible, by way 
of illustration, to calculate the flux referred to a room tempera- 
ture, regaixiless of o.ctual average temperature, i,e, to multi- 
ply "n" by 2.>. ,10-^cm/sec. The result is 

niox 

^ max 2.10^^cm“^sec"^ 

In practice it is rather difficult to shut down instsintaneously 
a pulsed reactor with large super criticality i,e, for the time 
smaller than If the reactor is shut down due to the negative 

linear thermal effect the reactivity decrease is proportional to 
th*- accumulated heat; drops four times. Besides, the real super- 
cr aility is always smaller than £«, -1, 

he neutron density can be inci^eased beyond 10^^cm“^ only in 
plasma-type reactors, i,e, by heating the core to 10^ deg, and above. 
However, the reactors of such types differ in principle from those 
which have do far been mastered by the reactor engineering* 

It may be mentioned that 6 /E^ is the ratio of the fissioned 
nit 
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nuclei t* the full number of the nuclei in the core 

the number of fisalon events per volume unit;. Consequently <S/E^<1. 
The ultimate limit, which is not achieved oven in the atom bomb 

B ~ 1025c.-5 

3. Heactors with tUel Kotion Belative to the Cozre 

in these reactors the fuel, thoroughly mixed with liquid or 
gaseous coolant, moves through the core at a high speed (homogeneous 
reactors). These reactors also include cyclopiles (Ref, 5) in which 
the zone with a high neutron flux moves along the long reactor 
body. Since it does not matter what is at rest - the geometrical 
region where fission is going on (core) or the mixture of the fuel 
with the coolant, this type of reactor actually differs only in 
the fact that the coolant may be in any aggregation state and, par- 
ticularly, in a solid state. 

Let us assume that a wave of the neutron flux and of the 
heating-up travels along the reactor with a velocity U, The absor- 
bers which were used to make the reactor subcritical are removed 
from the cold reactor body at the forward front of the wave. Absor- 
bers may be inserted into the hot core at the rear front of the 
wave to shut down the chain reaction ao as to avoid excessive hea- 
ting, If the negative thermal effect is sufficiently large there is 
no need to use rear absorbers. 

The homogeneous reactors are not usually provided with front 
or rear absorbers. Instead, the cross-section of the stream sharply 
narrowed, as it leaves the core, the neutron leakage increases enrt 
the reaction is shut down over a short section. 

For the rough description of the process given below the 
difference between the insertion of an absorber and narrowing of 
the cross-section is inessential. 

It would seem at a first glance that in this type of reactor 
the neutron density can be increased limitless even with a given 
heating limit by merely increasing the velocity U of the motion of 
the boundary which divides the region of the cold subcritical reactor 
from the supercritical region of the core, 

nk 


4 - 
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Of course, this involves engineering problems - the tested ma- 
terials should steadily keep up with the travelling neutron field, 
so •that they were subjected to irradiation during the time while the 
wave travels round the reactor body i if tho core is stationary the 
mixture of the fuel and coolant must be circulated at a high velo- 
city), Besides to preserve the duration of irradiation at the higher 
velocities the reactor length must be increased respectively (if 
the reactor is to cool dovm while the wave travells roimd the re- 
actor the capacity of the cooling system must grow with the velo- 
city U). However, these problems can be solved in principle. 

It appars, however, that as the velocity U increases the 
neutron density first grows in proportion to U but then tends to 
a certain limit. To illustrate this let uf compose a simplified 
unidimensional single-group equation which describes the distribu- 
tion of the neutron flux ^ along the corej 

']r' ( 7 ) 

Here - the difference between the material parameter 

and the radial geometric parameter V and D - the velosity and 
diffusion constant for the fission-producing neutrons. 

In the case of a travelling wave the time derivative ( ^ is 
the function x - ut) can be expressed throuf^ the coordinate 
derivative , 

at 

We obtain the equation 

li- 


^ 6 


d<f) 

dK 


( 8 ) 




( 9 ) 


U-x^ irj!> 

•% 

Outside the core dl? x is a hi^ negative value. Inside the 
core and YD are, generally speaking, variallee as the 

heating affects the neutron balance and the operating spectrum. 
This fact is very important when n^,^ has to be calculated for 
some specific reactor. To find the maximum value n^ ^ ^ they can be 
substituted with some average values. 

Let us first assume that in the reactor portions blackened 


by the absorbers the flux is zero, i.e. 

37<i 

- 5 - 
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Then the 
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critical len^h of the core (.the dlotance between the front and roar 
absortera ) is equal to 


where 


a 




n 


do; 


Oil) 


At 


y =1, obtain a^= co . The corresponding velocity U is yet small 


- r ^ - ^ w 

as compared with the neutron velocity ns D 36^ «• 1, At Y'' 1 ^q.(3) 
(9) has no solution that would become zox’o at two values of X. 

This case will be discussed separately. 

The distribution of the flux and the neutron density within 
the reactor is described by the formula 

^ ( 12 ) 

(X = 0 - rear boundary, X = a^ - front boundary of the reaction zone.) 

At small values of Y the cuirve is symmetrical, but when 

is close to 1 the maximum is displaced towards the rear front of 

the wave. The non-uniformity factor of the neutron flux distribution 

within the core is equal to 




rnojt 

n ev. 


T 


e- 

/ z ^~ ~ IT 


The integral neutron 
reactor can be represented in the foxin 


(15) 

ftr stme point of the 


X 


I n oLi ~ rid V. 

where Cl% / U/ - the time taken by the neutron wave to pass the 

point. Substituting with we obtain 


J 


“** = W,£f 


where 




(14) 


(15) 


^eff - UAL>, 

The fomnila (14) is similar to (2), Therefore, the equation for 
“max obtained by substituting in Eq.(5) with - 

the effective width of the curve n (b): 


3711 


(16) 
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Substituting ^ ^ 1*1(1. (16) with the cerresptndlng oxi)ro8si«n 
given by Eq.(15), substituting U in accordance with Kq,(lli and 
using Eq,(10) and (IJ) we obtain 


n 


ntnx 


(i/f 


i/f 


dil ff inji 


W(y) 




where 


yfy) ^ y 


i^C- 




■ T 




(17) 


(18) 


to equal jr at 
is below iznity and becomes 


Function 'X' monotcnically increaeea with 'y 
y" =1, Factor ■ 

unity only if coincides with the material parameter (very 

great transverse dimensions) and the migration area is fully 
described by D/£ , 

Let us now take up the case when X > 1. We will write the 

solution in the form 

It does not become zero at X=0 or X=a^. However, if S 
sufficiently small and a^ sufficiently great, n will be small 
at both boundaries of the core. This is in accordance with the 
actual condition that in the blackened areas n 0. The neutron 
density here depends on power of the available neutron sources 
(spontaneous fission, dxtemal irradiation, at the rear front- 
delayed neutrons) and on the degree of subcriticality. Naturally, 
it is very small compared with n^^^, especially at the forward froit, 
Asymptotically Fiuiction (19) has an exponential form 


(19) 


is 


exp 




Therefore, if at X = the flux is to be by several orders 
smaller than max ( ) - V 7 ^- / ) should 

he equal to at least a two-digit number. 

The non-uniformity factor is now described by the formula 


( 20 ) 


nil 


- 7 - 


Approved For Release 2009/08/26 : CIA-RDP88-00904R000100110021-6 





Approved For Release 2009/08/26 : CIA-RDP88-00904R000100110021-6 


Foi’ we obtain the same formula (17) in which, however 

'f'M j, ^ / 

At y =1, 'f' remains equal to , As Y rises the fun- 
ction V' omnotonically increases tending to unity at Y - > ' 

Therefore, the maximum neutron density (U — -y *-Kwl ) 

does not differ from the maacimum neutron density for pulsed re- 
actors, Of coui’se, instantaneous shut down of the reaction by the 
rear absorbers ( ^ =0) is 6m idealization, as is the instantaneous 

shut down of a pulsed reactor. 

As is the case with pulsed reactor hjjjgj decreases if negative 
thermal effect (decrease of during heat— up) is used instead 

of forced shut down. 

4. Reactors with Heat Transfer from Fuel to Coolant 

In reactors of such type high neutron density can be obtained 
only with small transverse dimensions of the fuel elements, minimum 
volume of the structural elements and small transverse dimensions 
of channels with the coolant, therefore, physically this reactor 
is almost homogeneous type, let us represent this reactor schema- 
tically in the form of a miarture of nuclei of two kinds. Fuel nuclei 
with atomic weight are as a rule at rest, macroscopic scatte- 
ring cross-section in the fuel elements is equal to and 

their volume part in the reactor is The average cosine of the 
scattering angle in the centre-of-mass system ±3 For the 

coolant whose nuclei move at an average velocity of corres- 

ponding characteristics are designated by index 2. Then the values 




O - O - O . 

^ ■ ,% c- — > 'i.-y n- — — 3 - (22) 

will be the probabilities of scattering on the fuel and coolant 
nuclei. 

The velocity of a newly formed neutron may be directed to 
any side with equal probability. Therefore the average value of 
the vector is equal to zero. Assming that ^ are indepen- 

dent from the neutron energy it is quite easy to calculate the 
average value of the vector Vj^ - the velocity of the neutron after 
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/ r .■f'j 


- H 


K collialon events (See the; A|)|)ondlx): 

/' /-.' A) , /] // ' AO '^.v) / , 

f- ..'Ll — — Lj'LL ' - 25 ) 

I > /f, I * t. • 

The time taken by the several initial collisions of the neutron 
is usually short as compared, to "the neutron life. Thorefoi'c, prac- 
tically during its entire life the neutron in general moves in the 
direction of the coolant flow (i.e. in the negative direction) at 
a velocity the absolute value of which es equal to 

V IV'-* 

^ n r. ... n \ i\ f I .. M. fi ,\ r Oil \ 


i- ^ ( 2 A) 

/ -f Jht I -f- 4! i 

k& U 2 is small as compared t# V, the variation of the absolute velo- 
city of a neutron with time practically is similar to that with 
slowing down in a stationary medium. Consequently motion of a part 
of the medium can be well evaluated by evaluating the averaged mo- 
tion of the neutrons at a velocity U^y» ^ additional term will 
appear in the balance equation of a stationary reactorj which will 
be written in the right-hand peart t 

(a similar equation is treated in Ref ,6), 

Equation i25) is obtained from (9) ty substituting U with 
Therefore it is desirable to redetermine the parameter If : 

However, the formula (15) for still contains U 2 . As a 

result an additional factor-ratio Ug/Ug^y appears in the formula (17). 
SlMiltaneously, substituted with SgQ^^q^il.e, heat 

accumulated by the coolant has to be distributea throxighout the 
entire c«re# We ©btain 


- 9 - 
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II fly,--) 


C ^ iHbaorbtioii cross-sect ion in the core and not in the 

coolant). To enable the reaction to be maintained at 1 the 

coolant entering tlie reactor must cany at least small quantities 
of neutrons. Othoi^wise, even with y -1 the reaction can not be 
maintained in the reactor with a finite height, of the core. 

To sum up, the maximum neutron density in the reactor with sta- 
tionary fuel and moving coolant differs from the maximum density in 
the homogeneous or pulsed reactor by an additional factor 
In accordance with (22) and (2^0 this factor is equal to 

r, r /Is, .)! i- 

-f— -/)/ — u.’ti _ 

lioA/ L JJ l-h/it.fli. 

itFx 




If S;j0 the factor is equal to unity, as would be expected (bomo- 
geneous reactor). It is also equal to unity when ^ 

are equal to ^ respectively. In the general oase^it 

may differ somewhat from unity, but it is very doubtful that any 
substantial gain can be obtained from the stationary condition of 
the fuel (it is not worthwhile to increase at the 

expense of low density of the coolant since also decreases). 

5. Reactors with Opposing Streams 
The slower rate of growth of n when the coolant velocity 
approaches 2TDSC.^ is due to the increase in the length of the region 
where the chain reacton is going on; since with the lengh of 

the fission reaction zone rises in proportion to U the heat liberation 

per volume unit of the core no longer grows, i.e, the neutron densi- 
ty n reaches its maximum limit determined by the formula (^), To 
overcome this limit the amount of neutrons carried away from the 
reactor by the coolanb must be decreased. As was kindly advised to 
the author by N.N.Ponomarev-Stepnoy this can be achieved by opliuting 
the coolant flow into stiteams and directing the adjacent streams so 
that they face each other. I«t us find the required degree of the 
coolant flow splitting. The reactor to be assumed homogeneous. The 
streams will have rectangular sectiwa with sides q, and q„ times 
374 ” ^ 

- 10 - 
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smaller than the corre spending dimonaienB ef the reactor a_ 

If the velocities of the streams +U^, then, using the Fourier 
series for a "chessboard” wo obtain: 

< - //„! r fz] 

This formula must be substituted into the equation 

(30) 

which is true if ^j/lj and great as compared to D, 

Let us seek the solution of (50) in the form of a series by 
the inherent functions of "an unperturbed" operator 


having the form 




The corresponding inherent values are equal to 


^ ^ - Sse^- 


where 


S3& 


l-af ^ a^) 


is the variation of the critical material parameter caused by "per- 
turbation" - motion of the coolant (the perturbation operator is 
equal to ). 

For expansion factors 4> by n,l,k obtain, in the 
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usual way a sot of oquatioiia in tho rom: 
^ +21 />„ ■ n, o C 

nnAiC U. 


*n},P,u 


0 


( 35 ) 


The matrix elomont contained in this equation is determined by the 
relation 

S/J, 


r ff'KU, % % ( 56 ) 


m n, 

where 


/\ (h) = 


I IfAis an odd number 

if h is an even number (57) 

The equations become simpler at great q^.q^, when approximately 

T_ Uo 

^nJtK •/ 


^ __ ctx ]r^ «o CO j 

■be easily seen that if ' -> 


^n,t,k. 


It ca^ be easily seen that if values n. I ( 2 s+ 1 )a i / id,. .. 

are .f the .ider .f the unit^ then ^ ^ 

jrry£ ^ _ 

» ".l-k are ef tta .laer ef the uni'ly, then ' 

- 12 - 
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If n=l=k=1, only thooo torms aro OQoontlal in the right-hond 
part of 08) in which m., j C2Q^-1)qy-P j , |( 2 t+ 1 )qj^-U | is of the 
order of the unity. 

Bearing in mind bhesa considerations the following formula can 
be easily obtained to change dC^ compensating for the entrainment 
and carry-over of the neutrons by the flowi 

f - rg y _ . 




where 


str “ 


S,ir3> £sh 








is the geometric parameter of one stream. 

The factor in square brackets in O^) does not practically 
differ from unity so that the effect of the neutron entrainment 
on the critical size is insignificant if "y < 1 , 

Disregarding the factor which is close to unity and consi- 
dering (54) we write (59) as follows: 




where 




i-h 


~y' z. ^ 

0 ^ 3i/a3S 




Cl!% 

The value q characterizea the degree 


of splitting the flow 


into streams. For each value of the parameter ( ^ 0 / 4 ^) there is 
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an optimum value of at which u,^ ia closoot to tho mini- 

mum value of /X, i.G. the neutron doriuity ie maximum. 'I'o en- 
aui'e that the difference botwcon find 'J' /<£. ia not yi'cat 

the parameter C Aq) exceed unity, Wlienco the 

X ) 

required degree oi splitting : 

For example, to raise the neutron density by one order higher ns to 
its maximum value in a pulsed reactor must be approximately 

15» Q qi approximately 4, 

6. Conclusion 

It is not an easy problem to obtain neutron densities dose 
to the maximum values determined by the formulas (4) and (27). 

With fast— neutron reactors this becomes downright impossible 
since there are no methods by which the reaction can be shut down 
during the time determinea by the formula (1). Even scattering 

of the materials due to their thermal expansion or evaporation does 
not proceed quickly enough (naturally, at ^ <C 1 ev). The theore- 
tical limit can be reached more easily in a pulsed resonance reactor. 

With thermal-neutron reactors the theoretical limit can be 
closely approached not only in a pulsed reactor but also in a homo- 
geneous reactor. The required velocities of relative motion of the 
core and the reactor body come close to 10^ m/sec (at ^ ^ 1). 

For reactors in which the fuel is separated from the coolant such 
velocities are less feasible than for homogeneous reactors as the 
fuel elements must withstand the coolant flow but their thickness 
must be reduced to the minimum to avoid excessive internal overheating 
Pu 37 ther inci^ease of the velocity (U > 2VD X ) will be of no avail 
unless opposing streams of the coolant are used. In a cyclopile 
this complication of design is very important ♦ So, inside the re- 
actor with a non-ionized coolant the neutron density of order 
lO^'^cm”^ can be obtained only in case the neutron spectrum is thermal. 


^Strictly speaking with ^ y4q ^ 1 formula (59) cannot be 
applied. 


in 


- 14 - = 
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However in the reflector a higher neutron density than 10 cm 
is avllable hy forcing the resonant reactor. 

With resonant spectrum y = 1» wtien U is measured by kilo- 
meters in second. Therefore resonant neutron carryover is not 
important, the length of the core of a resonant reactor is inde- 
pendent on U and thermal stress is limited only by a permissible 
velocity of the fuel and the core motion. The neutron density in 
the reflector, static to the core, does not directly depend on fuel 
concentration in the core, l»e, hardness of spectrum. This density 
may be 1-2 orders more the neutron density in the core of the re- 
sonant reactor and the limit established by formula (4) is not 
valid for it. 

But one cannot count on a high gain compared to the thermal 
reactor. In the reflector heat is liberated which is to be removed. 

So thermal neutron carryover by a coolant is also limiting 
the available lux. 

For many purposes, merely obtaining high neutron density 
is not enough; this density must be retained for a comparatively 
long period of time. Obviously, in this respect the cyclopile or 
homogeneous reactor are superior to the pulsed reactors. It also 
quite clear that production of a high neutron flux to be maintained 
for long periods of time with minimum expenditures places specific 
requirements on the selection of the core materials^ its size and 
construction. 

7. Appendix 

Assume that and ?' are the neutron velocities in the 
c c 

centre-of-niass system tefore and after collision with a nucleus 

whose velocity and mass number in a laboratory coordinate system 

are u and A, respectively. As 7^ = it can be ass^xmed that 

results from the effect produced on V by some random turn ope- 

c 

rat or i.e. 

The average cosine of the turn angle is ^ i i#e* 


C'l) 

(. 2 ) 
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As the average value of the vector 0 V 

c 

is the vector parallel to V , from (2J it follows 

) a./. Vc 

^ ^ ^ (3) 

Take notice, ^ is a linear operator. 

The neutron velocities in the laboratory system before and after 
collision (V and V* ) are related as follows 

iri tiM 

l-f-/! ^ l-^/l l-t-fi ^ jirfi ^ 

If there are two kinds of nuclei and the probabilities of 
collision with them are equal to P^and Pg then 

V' = 17 + 0 (5) 

where L and 0 with probabilities are respectively equal to 

ITJi ( 6 ) 

'i-hJh 


From the formula C5) it follow* that by averaging (5) with res- 
pect to collision the operator S can be substituted with ^ 


so that 


where 


V ss L 7 + 0 
av av av 


^av “ ^1^av ^2^av* ®av ~ ^1®1av ^2®2av 

oM and are obtained frwm (6) and (7) "by substitu- 
ting with, i-*®* 






“av irkr^^ 


( 10 ) 


How, let us take a neutron with a randomly oriented velocity 7^ 

after K collisions. Its velocity 7. can be obtained by applying K 
374 ^ 

-16 
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times the formula The average value is obtained by 

subsequent averaging for all collisions and averaging for initial 
direction V^. However, since the collisions are statistically 
independent (we take that and y'Vi ^ are independent on energy) 
the same result can be reached by applying K times the formula 
(8), assuming that „„ = 0. We obtain: 

O civ 



in 


-17 ~ 
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